Stable Propagation of a Burst Through a One-Dimensional Homogeneous Excitatory 

Chain Model of Songbird Nucleus HVC 
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We demonstrate numerically that a brief burst consisting of two to six spikes can propagate in a 
stable manner through a one- dimensional homogeneous feedforward chain of non-bursting neurons 
with excitatory synaptic connections. Our results are obtained for two kinds of neuronal models, 
leaky integrate-and-fire (LIF) neurons and Hodgkin-Huxley (HH) neurons with five conductances. 
Over a range of parameters such as the maximum synaptic conductance, both kinds of chains are 
found to have multiple attractors of propagating bursts, with each attractor being distinguished 
by the number of spikes and total duration of the propagating burst. These results make plausi- 
ble the hypothesis that sparse precisely-timed sequential bursts observed in projection neurons of 
nucleus HVC of a singing zebra finch are intrinsic and causally related. 
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I. INTRODUCTION 



A question of great interest to neurobiology is how an- 
imals learn to generate temporal patterns of muscle acti- 
vation. An example that has been much studied in recent 
years because of its relevance to human speech Q and be- 
cause of the rich variety of possible experiments is how 
songbirds learn to sing a song by auditory-guided vocal 
feedback . A young male bird first memorizes the song 
of an adult male of the same species. Then over many 
months, over many iterations (more than 50,000 itera- 
tions for a zebra finch Q), and by just listening to his 
own vocalizations, a young male learns to vary the activa- 
tion of its respiratory and syringeal muscles until his song 
is able to match accurately the original memorized song. 
This is an impressive feat given that an adult song might 
last several seconds, that there is auditory structure that 
lasts less than 10 ms, that many syringeal and respiratory 
muscles need to be coordinated, and that some species 
of songbirds are able to learn and sing many different 
songs. Nearly all details of this process remain poorly 
understood, in particular how the young male memorizes 
a complex song, and how auditory feedback is used to ad- 
just the pattern of muscle activation until the songbird 
is able to reproduce accurately the memorized song. 

Studies of songbirds have shown that certain anatomi- 
cally and physiologically distinct brain regions called nu- 
clei are associated with the recognition, learning, and 
production of song (see Fig.Hfr. A recent experiment by 
Hahnloser et al [j] recorded extracellular action poten- 
tials (spikes) from neurons in awake singing male zebra 
finches and found that the neurons in the nucleus HVC 
that project to the robust nucleus of the arcopallium (ab- 
breviated as RA p|) have the remarkable properties of 
firing sparsely and precisely during singing. (In the fol- 
lowing, HVC neurons that project to RA will be abbrevi- 
ated as HVCra neurons.) Typically, each HVCra neu- 
ron fires a brief burst of about 3-4 spikes once per song 
motif with each burst lasting about 6 ms. (A motif is a 




FIG. 1: Schematic anatomical diagram of a sagittal cross sec- 
tion of one hemisphere of a zebra finch brain, showing the ma- 
jor nuclei associated with the learning and production of song. 
The nucleus HVC plays a central role since it receives auditory 
and other input, it drives the premotor nucleus RA Q, and 
it sends information to a parallel anterior forebrain pathway 
(nuclei X, DLM, and LMAN) that play a role in the learning 
and maintenance of a song. There are two HVC nuclei, one 
on the left and one on the right side of the songbird's brain. 



cluster of distinct auditory syllables that is repeated as 
a single pattern and that lasts an average of 0.6 s for ze- 
bra finch songs.) Measurements during successive motifs 
from a given adult male bird show that each burst from a 
particular HVCra neuron is aligned with certain acous- 
tic features of the motif to a precision of about 0.7 ms Q. 
These bursts are important to understand since they are 
believed to provide the temporal framework for organiz- 
ing the syllables of a song jg] . 

The observation that each HVCra neuron bursts just 
once per motif and other data which show that there 
is a reproducible temporal ordering of HVCra neurons 
according to when they fire their single burst per mo- 
tif 3 , suggest the hypothesis that bursts propagate along 
a feedforward network of excitatory HVCra neurons [|| . 
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The idea is similar to a path of dominoes such that the 
falling of one domino causes the next domino in the path 
to topple. Here, a given HVCra neuron receives one or 
more bursts from other HVCra neurons that have fired 
recently during the motif. In response, the given neuron 
itself fires a burst that helps to activate neurons that fire 
later during the motif. That each HVCra neuron fires 
exactly once per motif implies that the connections are 
feedforward (although this assumption can be weakened 
if a recurrent burst arrives while some earlier neuron is 
still in a refractory state 0). The bursts continue to 
propagate until the chain ends and the motif stops, or 
perhaps there is a branch point to a different chain whose 
dynamics generate a different motif. 

While there has been much previous theoretical re- 
search concerning when spikes can pro pagate along exci- 
tatory networks of different kinds [lU (13111 El El 0, 
the songbird experiments on HVC of Hahnloscr et al dur- 
ing singing Q and related experiments by Mooney of 
bursts in HVC projection neurons during audition Ell 
motivate asking new theoretical questions about the 
propagation of a burst through various kinds of networks. 

In this paper, we show that a minimal feedforward 
architecture — a one-dimensional (Id), homogeneous, ex- 
citatory chain of non-bursting neurons — can support the 
stable propagation of a brief high-frequency burst. By 
studying two kinds of chains, one with leaky intcgratc- 
and-fire neurons (abbreviated as LIF neurons) and one 
with more realistic single-compartment conductance- 
based neurons (abbreviated as HH neurons for Hodgkin- 
Huxley like), we show that a burst can propagate in a 
stable manner that does not require a careful choice of 
neuronal model nor a careful tuning of model parame- 
ters. Depending on the values of parameters such as the 
maximum synaptic conductance, we find that each kind 
of Id chain (LIF or HH) has multiple attractors that dif- 
fer in the number of spikes and in the total width of the 
propagating burst. 

Our calculations are the first to demonstrate that 
a brief high-frequency burst similar to those observed 
in HVC can propagate in a stable manner along a sim- 
ple excitatory chain of neurons. (A similar result was 
announced independently by Jin and collaborators Q, 
who studied a more complicated model. We discuss 
this preprint briefly in Section llVI below.) Our results 
make plausible one of the simplest explanations of the 
Hahnloser et al result, namely that the observed bursts 
of HVCra neurons during singing are intrinsic to HVC 
in that external input from other brain regions is not 
needed to generate the bursts (except the first burst), 
and in that the bursts are causally related such that one 
burst initiates the generation of the next burst. Our re- 
sults also make some predictions such as the existence of 
multistability (different kinds of propagating bursts can 
occur in the same chain depending on how the chain is 
activated), and that transitions from one kind of prop- 
agating burst to another kind (differing in the number 
of spikes and burst width) can occur as parameters are 



varied. New experimental studies of HVC, especially us- 
ing optical methods with high time resolution that 
can examine the spatiotemporal activity of many neurons 
simultaneously, may help to confirm these ideas. 

The rest of this paper is organized as follows. SectionlTTI 
discusses some details of the two classes of mathematical 
neurons that we use in the Id chains, and how the excita- 
tory synaptic currents are modeled. Section IlIII presents 
numerical results mainly for a Id chain of HH neurons (as 
opposed to LIF neurons), especially the existence of dif- 
ferent attractors corresponding to different kinds of prop- 
agating bursts. The paper concludes in Section llVl where 
key results are summarized. We also discuss there how 
our results are related to alternative hypotheses such as 
that the bursting arises from propagation through synfirc 
chains |8| or by a central pattern generator (abbreviated 
asCPG). 

In a second paper [TtI ITsf , we will discuss how noise 
and network heterogeneities affect the propagation of 
bursts through a synfire chain. This second paper shows 
that our key conclusions still hold: a brief high-frequency 
burst can still propagate in a stable manner for a range 
of noise strengths and for various amounts of heterogene- 
ity. However, we also find that single spikes, even if syn- 
chronized across a synfire pool, do not propagate in a 
stable way for most of the parameters studied, which 
suggests that short bursts of several spikes are important 
for achieving robust propagation through realistic synfirc 
chains. 



II. METHODS 

In this section, we discuss some experimental data 
that justifies the study of a one-dimensional chain of 
excitatory neurons, and we also point out some of the 
experimental details that are ignored in our model. 
We also discuss some details of the integrate-and-firc 
and conductance-based neuronal models with excitatory 
synapses that we use in the Id chains. Appendix ^ pro- 
vides further details of equations and parameter values. 



A. A Homogeneous Id Chain of Excitatory 
Neurons 

The physiological properties of HVC neurons and how 
they are interconnected within HVC are poorly under- 
stood [H EM IMEl so that it is not possible at this 
time to develop a quantitative model of the HVC micro- 
circuitry, say at the level of hippocampus models [23j . 
Researchers have shown |0, [^J that there are at least 
three main classes of neurons: excitatory neurons that 
project to nucleus RA, excitatory neurons that project 
to area X in the anterior forebrain pathway, and in- 
hibitory interneurons [Tol l2l| that connect only to other 
neurons within HVC. Recent paired-electrode recordings 
by Mooney and Prather [20j have shown that each type 
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of HVC neuron makes local connections with the other 
two types of HVC neurons but details of the connections 
such as the number, kinds, and strengths of synapses are 
incompletely known. 

Of special importance for this paper is the experimen- 
tal observation that HVCra neurons synapse with other 
HVCra neurons |2jj. Thus it is possible for a feedfor- 
ward network of excitatory neurons to exist in HVC, al- 
though we emphasize that there is no evidence presently 
for such a network. Further, there are of order 40,000 
HVCra neurons in HVC p , which should be sufficiently 
many to create a chain whose dynamics spans a mo- 
tif or even several motifs. (Since the observed bursts 
in HVCra neurons last 6 ms Q, a chain of about one 
hundred pools of neurons would suffice to span a motif 
of 0.6 s, provided that the bursts do not overlap in time.) 

To determine in principle whether the propagation 
of a burst can explain the experimental data 0, Il5j . 
we study one of the simplest possible feedforward net- 
works, namely a one-dimensional homogeneous feedfor- 
ward chain of identical neurons such that each neuron 
connects via a single identical excitatory synapse with 
the next neuron of the chain. The assumption of homo- 
geneity is not realistic biologically but reduces the pa- 
rameter space of the calculations and increases our abil- 
ity to understand how the existence and properties of a 
propagating burst depend on parameters. 

Our Id model leaves out many experimental details (we 
discuss this further in Section ITTT1 below), of which two 
are especially important, the need for multiple pathways 
and the presence of inhibitory neurons. First, as pointed 
out by Abeles and others, transmission of information 
along a one-dimensional chain is not robust since damage 
to any neuron in the chain can alter or stop the trans- 
mission of information. (That HVC is robust is demon- 
strated by the fact that the adult song of a zebra finch 
changes little over the life of the bird, even though there 
is ongoing neuro gene sis and hence a steady turnover of 
HVCra neurons |24 l25f.1 If some kind of propagation 
occurs in HVC, there must be parallel pathways. 

A widely studied example of a feedforward neuronal 
architecture with multiple pathways is a synfire chain 0, 
U |2(|, see Fig. |2Ja) . Although there are variations of 
this architecture, a synfire chain is usually described as a 
strictly feedforward chain of pools of neurons such that 
the neurons in one pool project only to neurons of the 
next pool. Reliable transmission occurs even in the pres- 
ence of noise or of network heterogeneities because of 
the multiple paths, and also because each neuron re- 
ceives spikes from many neurons of the previous pool and 
the spikes spontaneously become highly synchronized as 
propagation continues ja • The arrival of synchronized 
spikes at the many synapses of a given synfire neuron 
causes a large postsynaptic current that triggers new 
spikes with a probability that can approach certainty. 

The idea that the sparse brief bursts observed in 
HVCra neurons are causally connected by sequential fir- 
ing through some kind of chain has been proposed by 
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FIG. 2: (a) Schematic diagram of a feedforward synfire chain 
consisting of successive pools of neurons (vertical column of 
circles). Each pool contains excitatory HVCra neurons that 
fire approximately synchronously to activate neurons in the 
next pool. The horizontal hollow arrow on the left represents 
input that can initiate the synfire chain. The hollow arrows at 
the top of each pool indicate efferents that could convey infor- 
mation from a given pool to other brain areas, (b) Schematic 
diagram of the one-dimensional homogeneous feedforward ex- 
citatory chain that we study in this paper. In the biologically 
unrealistic but theoretically convenient case that all neurons 
in the synfire chain of (a) are identical, that all neurons in 
one pool connect to neurons in the next pool in the same way 
with identical synapses, that axonal delays between pools can 
be ignored, and there is no noise, the ld-chain will have dy- 
namics identical to a synfire chain whose spikes have become 
fully synchronized. 



experimentalists [ a. [271 and assumed in several theoreti- 
cal studies j2^, I29L l30| . In a recent experiment [2(j , cal- 
cium imaging and timing studies of intracellular record- 
ings have suggested the existence of synfire chains in the 
visual cortices of mice and cats but no similar imaging 
study has yet been carried out for songbirds. 

The dynamics of our onc-dimensional homogeneous 
chain is equivalent to the dynamics of a homogeneous 
noiseless synfire chain (homogeneity here means equal 
numbers of identical neurons per pool all connected in 
identical ways from one pool to the next) such that all 
spikes within a given pool have become perfectly synchro- 
nized. This is an important observation since the synapse 
that connects one neuron to the next in the Id chain of 
Fig- IHb) then has to be strong to correspond to the many 
synapses of a neuron in the full synfire chain. In general, 
noise, heterogeneities, delay times, and imperfect syn- 
chronization of spikes cause the dynamics of a synfire 
chain to differ from that of our Id chain so that simu- 
lations of real synfire chains are needed to understand 
these effects on the propagation of a burst @, El El • 

A second important experimental detail that we leave 
out of our model is the presence of the inhibitory HVC in- 
terneurons. Our justification for this assumption is sim- 
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ply that the function of these interneurons is not known. 
Proposed mathematical models for sequence generation 
and sequence recognition incorporate inhibitory neurons 
in various ways [3^, |3^ but there is no support 

yet for these models in the context of the songbird sys- 
tem. We argue that ignoring the inhibitory neurons is 
a useful first step toward understanding the HVC mi- 
crocircuitry since the calculations discussed below show 
that a minimal idealized Id chain of identical excitatory 
neurons already suffices to produce a stable propagating 
burst, in which case other effects such as noise, inho- 
mogeneities, and interneurons could be included pertur- 
batively. Further, some experiments suggest that ignor- 
ing the interneurons may not be too severe an assump- 
tion. HVC interneurons tend to fire rapidly, chronically, 
and in approximate synchrony throughout a motif [6j, 
and rapidly and chronically during audition of a bird's 
own song [15], so that the interneuron spikes might not 
be closely correlated with the brief sparse bursts of the 
HVCra neurons. In this case, the main purpose of the 
interneurons might to prevent runaway excitation of the 
projection neurons. 

We conclude this section with the technical observation 
that our assumption that the Id chain is homogeneous 
and that there is no delay along axons (the axons along 
the chain have effectively zero length) allow the Id chain 
to be simulated by integrating numerically a single neu- 
ron with a single afferent synapse. During the integra- 
tion, the output of the neuron can be stored in memory 
and then used as input to the same neuron which then 
represents the next neuron in a chain. A related idea was 
used by Reyes [TTI |. who used a single real hippocampal 
neuron in a slice together with the dynamic clamp tech- 
nique to simulate a synfire chain consisting of many pools 
of identical neurons. 



B. A Leaky Integrate-and-Fire (LIF) Model of an 
HVCra Neuron 

We studied propagation of bursts along Id chains con- 
sisting of one of two types of neurons, leaky integrate- 
and-fire (LIF) neurons that are described in this section, 
and more realistic conductance-based Hodgkin-Huxley 
(HH) neurons that are described in the next section. Nei- 
ther neuronal model bursts intrinsically when subject to a 
direct current (DC) external stimulus. Provided that an 
initial burst is applied to the first neuron of the Id chain, 
we find that successive neurons are capable of generat- 
ing an identical output burst when stimulated itself by a 
burst. 

With respect to conductance-based models, LIF mod- 
els (34) are attractive since they have fewer parameters, 
lead to efficient numerical simulations, and are sometimes 
amenable to analytical methods. The fcth LIF neuron of 
the Id chain (fc > 1) was modeled as a first-order non- 



autonomous ordinary differential equation 

T m — = v ~v k + R(l k (t) + I k (t)), (1) 

where the parameter r rn is the capacitive i?C-time scale 
of the neuronal membrane, the variable v k (t) is the po- 
tential difference across the membrane of the fcth neuron 
in the chain, the external current I k {t) of the fcth neu- 
ron is a specified function (in this paper applied only 
to the first neuron to initiate activity in the chain), the 
synaptic current Ig(t) arises from spikes in the previous 
(k — l)st neuron of the chain (see Eqs. and be- 
low), the parameter vq is the resting potential toward 
which the potential v asymptotes in the absence of ex- 
ternal and synaptic currents (I e — Is — 0), and the 
parameter R is the total resistance of the membrane. 

Eq. Q was supplemented with the usual spiking rule 
that, whenever the potential v k (t) is increasing and 
crosses a specified threshold value ^thresh > ^o, the neu- 
ron is assumed to have spiked and the potential v k is 
instantly and discontinuously decreased to a reset value 
Preset < v o- An absolute refractory period was included 
in this LIF model by freezing the potential v k to the 
value f r eset for a time interval Refract after a spike. We 
found that including a refractory interval of i ro fract = 
1 ms produced only minor qualitative changes to the re- 
sults obtained in the absence of a refractory period, for 
example a 1 ms refractory period expanded by a modest 
amount the basins of attraction for the different stable 
propagating bursts discussed in Section lTTTl fsee Fig. Elbe- 
low). Although HVCra neurons show some accommoda- 
tion 0| , we did not include this detail in our LIF model. 

The total postsynaptic current /| (i) of the fcth neuron 
in the chain was assumed to be a linear accumulation 

iIv^y, 1 ^-^ 1 ), (2) 

i 

of postsynaptic currents I s (t — t k ~ 1 ) associated with 
spikes that occurred in the previous (fc — l)st neuron of 
the chain at times t^ 1 < t. (Eq. J2J also assumes that 
there are no time delays arising from the time for a spike 
to propagate from one neuron to the next.) The time 
dependence of the synaptic conductivity associated with 
a single spike at time t k — was modeled as a double 
exponential [34|. leading to the following expression for 
the postsynaptic current per spike: 

I s (t,v)=nl (e~ t / T i -e-*A*). (3) 

The slow and fast time constants t\ and T2 had respec- 
tively the values 1.1 ms and 0.2 ms to match roughly 
the 4.0 ms onset-to-peak time constants observed in 
paired intracellular recordings of HVCra neurons pof . 
The value Jo = 0.3 nA was chosen to give a 1.0 mV peak 
excitatory post-synaptic potential (EPSP) which is con- 
sistent with experiment J20j . The number n, which varied 
from 1 to 32, indicated the synaptic strength in terms of 
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the number of synchronous spikes that a neuron would 
receive if placed in a uniform synfire chain with n neu- 
rons per pool and all neurons of one pool connecting to 
all neurons in the next pool (with synapses of the same 
form Eq. ©). The range 1 to 32 was found empirically 
to span the range of stable dynamics, see Fig. below. 

An homogeneous chain was then obtained by using for 
each LIF neuron the same ten parameters r m , vq, R, 
^thresh, Preset, Refract, n, I , n, and t 2 , and the same 
function Eq. © for the postsynaptic current per spike. 
A time constant r m = 15 ms was used to approximate 
the in vivo response time |l5j| . and a membrane resis- 
tance of R — 60 Mfl was chosen to match the impedance 
of HVCra neurons ^5|- Other parameter values used 
were v Q = -70 mV, v thlcsh = -55 mV, v rcsct = -75 mV, 
and tj-efract = 1.0 ms. 

The chain of LIF neurons represented by Eqs. dj)-© 
was integrated by using a forward-Euler method [35j with 
a constant time step of At = 0.01 ms. (Smaller time 
steps by factors of four did not lead to significant changes 
in the results.) The supplementary reset and refractory 
rules were applied at the end of each time step. For most 
runs, a zero external current I k = was assumed for 
each neuron of the chain except for the first neuron, for 
which a step function was used to stimulate a burst. The 
code was programmed and run using the computational 
mathematics program Matlab version 6.5 |36j |. 



C. A Single-Compartment Hodgkin-Huxley Model 
of an HVCra Neuron 

The second neuronal model that we used to study the 
propagation of a burst in a Id chain of excitatory neurons 
was a single-compartment model based on the Hodgkin- 
Huxley equations with five representative conductances. 
The evolution equation for the membrane potential v k (t) 
of the A:th neuron in the chain was 

£ 9i(t, v k ) {vi - v k ) + I k e (t) + I k s (t). (4) 

i—l 

The symbols have the following meanings: C m is the total 
membrane capacitance, gi(t,v) is the voltage-dependent 
conductance of the ith kind of membrane channel, Vi is 
the resting potential for the ith channel, and the cur- 
rents I k (t) and Ig(t) have the same meaning as in Eq. 
Appendix ^ gives the details of parameter values and 
other evolution equations related to the conductances gi. 

Although more realistic in terms of its time depen- 
dence, the HH model Eq. Q is not necessarily more sci- 
entifically appropriate than LIF models since the proper- 
ties of the various membrane conductances are only par- 
tially known for HVC neurons [3 13 IH E3, the spatial 
distribution of the channels in HVC neurons has not been 
determined (which would be needed to construct accurate 
multi-compartment models), and little is known about 
the type, strengths, and locations of synapses in HVC. 



We did find that a single-compartment HH neuron with 
five conductances is only able to match some features of 
the experimental HVC data, and it is not clear which 
experimental details are more important than others to 
include in the process of fitting a mathematical model 
to data. For these reasons, and because this paper is 
concerned whether in principle an HVC-like burst can 
propagate stably through an excitatory chain, the con- 
ductances of the HH model Eq. (0J were only loosely 
based on existing HVC data. 

Our starting point for choosing the conductances in 
Eq. was a recent paper by Prinz et al [3^, which 
showed that an eight-dimensional phase space obtained 
by varying the maximum conductances of eight channels 
(whose functional properties were obtained from lobster 
stomatogastric neuronal data) contained a great diversity 
of dynamical behavior. We chose the functional forms of 
the leakage, sodium, and potassium channels from this 
paper and added two other channels as suggested by 
current-clamp data of HVC neurons [33 : a fourth chan- 
nel was a low-threshold transient calcium current, and 
a fifth channel was a calcium-activated potassium cur- 
rent (denoted in the following by the symbol KCa). The 
calcium channel was chosen to activate around -50 mV 
and inactivate around -80 mV with time scales such that 
the channel could be activated transiently [3H l39| . The 
KCa conductance activates the potassium current after 
intracellular calcium concentration rises, with a func- 
tional form adopted from Yamada's model jijj. The 
maximum conductances of all five conductances were ad- 
justed iteratively by hand to produce a short spike width 
and to generate a spike adaptation with a high initial fir- 
ing rate as observed by Kubota 37] . A reasonably good 
fit to the Mooney in vivo data (see Fig. 1 of Ref. [Tj|) un- 
fortunately could not also match Kubota's in vitro spike 
profiles of Ref and vice versa. (We hope to obtain 
more complete fits to HVC data by adding more chan- 
nels and more compartments when more anatomical and 
physiological data become available.) 

The total synaptic current was also assumed to sat- 
isfy the linear relation Eq. J3J) but an alpha function [34[ 
was used instead of Eq. © to approximate the time- 
dependent probability of a channel opening. The exci- 
tatory postsynaptic current (EPSC) arising from a spike 
that occurs at time t = had the form 

I s (t;v) = g s C (J^\ e-t/r* (v s -v), (5) 

where the constant C with value e normalizes the max- 
imum value of the expression C(i/r s ) exp(— t/r s ) to 1. 
This normalization makes the parameter g s the maxi- 
mum conductance, which we varied over the range to 
0.55 nS. In Eq. (jSJ, the synaptic time constant r s was 
fixed with value 7 ms and the synaptic reversal volt- 
age v s was set to 120 mV to make the synapse excitatory. 
Eqs. (0J, @ and (JSJ) together with the evolution equa- 
tions for gate variables given in Appendix [S] were inte- 
grated with the Neuron simulation code version 5.6 |4l| . 
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with a constant time step of At = 0.1 ms. 



III. RESULTS AND DISCUSSION 

A. Results for a Id Homogeneous Chain of 
HH Neurons 

The calculations discussed below of one-dimensional 
homogeneous chains using the LIF or HH neurons of the 
previous section with excitatory synapses show that a 
brief high-frequency burst similar to that observed in 
HVCra neurons during singing 4] can propagate sta- 
bly over a range of parameter values. These results 
make more plausible the hypothesis that an excitatory 
chain of HVCra neurons generates the sparse precisely- 
aligned bursts observed experimentally [J. Our results 
show further that different asymptotic attractors can ex- 
ist for given parameter values so that hysteresis can oc- 
cur (different initial conditions can lead to different non- 
transient dynamics). The fact that the bursts exist as 
an attractor means that there is a transient time dur- 
ing which properties of the burst such as the number of 
spikes and burst width (alternatively, the average burst 
frequency) evolve until the final stable values are ob- 
tained. This transient time varies with the initial state 
used to start the chain and with the choice of parameters. 

Since these qualitative conclusions turned out to be 
similar for LIF and HH neuronal models, we report here 
results mainly for the HH models of Section III CI and 
mention briefly in the next section how the results differ 
for a chain of LIF neurons. Since the parameter spaces 
for LIF and HH neurons are high-dimensional (10 pa- 
rameters for LIF, about 25 parameters for the HH neu- 
rons, and these do not include the choice of the func- 
tional form I s (t, v) for the EPSC per spike), to establish 
our key results we varied only one parameter systemat- 
ically, namely the maximum synaptic conductance g s of 
Eq. JSJ for HH neurons, or equivalently the number of 
synchronous presynaptic spikes n of Eq. (J3J) for LIF neu- 
rons. 

Figure |3 shows how the same initial burst of five spikes 
propagates through a Id excitatory chain of HH neurons 
for several different values of the maximum synaptic con- 
ductance g s . The initial burst was created by using an 
external current ll{t) of Eq. J3J for the first neuron in 
the chain to inject a square pulse of current that caused 
five spikes to appear in rapid succession. For a synaptic 
coupling strength g s smaller than about 0.1 nS, Fig. [3^ 
shows that the initial burst rapidly dies out by the third 
neuron of the chain and all spikes disappear. Over a 
range of stronger couplings 0.2 < g s < 0.4 nS, the initial 
bursts evolves into a stable invariant propagating burst 
that can have one to five spikes (columns (b) and (c) 
of Fig. For stronger couplings g s > 0.4 nS, the initial 
burst is unstable and the number of spikes grows steadily 
without limit. However, other initial conditions can lead 
to stable bursts in this range, an example of hysteresis. 



(a) (b) (O (0) 

g =0 1 nS g =0.2nS g = 3 nS g = 45 nS 
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FIG. 3: Propagation of an initial 5-spike burst through an ho- 
mogeneous one-dimensional excitatory chain of HH neurons, 
for different maximum synaptic conductances g 3 of Eq. (y. 
Each column is corresponds to a fixed g s value, and suc- 
cessive rows of a given column show the membrane poten- 
tial v k (t) as a function of time for successive neurons in 
the chain (k = 1, 2, • ■ ■ ). (a) For a weak synaptic coupling 
strength g s = 0.1 nS, the initial burst decays by the third 
neuron, (b) For a stronger coupling g s — 0.2 nS, the initial 
bursts decays to an invariant single spike, (c) For g s — 0.3 nS, 
the initial burst evolves to a stable state with three unevenly 
spaced spikes, (d) For still stronger couplings g s > 0.45 nS, 
the initial burst can be unstable for some initial states and 
the number of spikes increases steadily. 



Fig. 01 provides a more global understanding of the 
transient dynamics and resulting attractors. In each 
panel, the vertical axis indicates the number of spikes 
observed in a burst while the horizontal axis indicates 
the position k of a neuron along the chain. For a weak 
synaptic coupling g s < 0.2 nS, all initial states decay to 
zero spikes (not shown in the figure). For a synaptic cou- 
pling of g s — 0.3 nS (Fig. 0k), initial bursts with three 
or more spikes decay within four neurons to a common 
final burst of three spikes. An initial burst with two 
spikes evolves slightly to an invariant burst with also two 
spikes, and a similar conclusion holds for an initial burst 
with a single spike. As the synaptic strength is increased 
to g s = 0.35 nS (Fig. HJa), the number of attractors in- 
creases so that stable bursts with one to five spikes are 
observed depending on the initial state of the first neu- 
ron. The transient time is still rather short, with the 
final number of spikes stabilizing in all cases within three 
neurons (k < 3). Finally for synaptic couplings stronger 
than about 0.4 nS (Fig.0Ji), stable and unstable states are 
observed depending on the initial state of the first neu- 
ron. The unstable bursts all grow at the same rate, with 
an extra spike appearing for each next neuron traversed. 

The transient dynamics of Fig.0]are presented in some- 
what finer detail in Fig. [S] which shows that a non- 
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FIG. 4: (Color onlin) The dependence of transient dynamics 
and resulting attractors of propagating bursts on the initial 
injecting current magnitude l\ (t) and on the synaptic conduc- 
tance g s . For weak couplings with g s < 0.2 nS, all initial states 
decay to zero (not shown here), (a) For g s = 0.3 nS, there 
are multiple basins of attraction corresponding to asymptotic 
propagating bursts with 1, 2, or 3 spikes. Transients decay 
rapidly by the third neuron of the chain, (b) Increasing the 
synaptic strength to g s — 0.35 nS increases the number of 
attractors, with final bursts containing 1 to 5 spikes, (c) 
For a larger synaptic conductance g s = 0.45 nS, initial condi- 
tions lead to stable or unstable states. An initial burst with 
6 spikes will evolve to a slightly different invariant burst with 
also 6 spikes (squares). 



transient propagating bursts is characterized by at least 
two parameters, their total width and the number of 
spikes. For example, the same beginning state B\ can 
evolve into three different end points depending on the 
synaptic coupling strength g s : point E\ with three spikes, 
point Ei with 4 spikes and with a slightly longer width 
of 45 ms, and an unstable state (triangles). 



Finally, Fig. [H] indicates the approximate basins of at- 
traction that exist for different initial current amplitudes 
(horizontal axis) and different maximum synaptic con- 
ductances (vertical axis). The plotted numbers indicate 
the number of spikes observed in the asymptotic state for 
the specified axis values. This plot shows again how all 
states decay away for sufficiently weak synaptic coupling, 
all states are unstable for sufficiently strong coupling, and 
that stable propagating bursts with 1 to 6 spikes can be 
found for intermediate coupling strengths. The basins of 
attraction turn out to be of comparable size except for 
the tiny basin corresponding to a 6-spike burst. 



FIG. 5: (Color online) Different synaptic conductances g a 
lead to different final propagating burst widths and spike 
numbers, which shows that the total burst duration is an- 
other important characteristic of the final attractor. For 
g s — 0.3 nS, the curve of block symbols evolves into the 
point El. For g s = 0.35 nS, the curve with diamond sym- 
bols ends in the point E2. For g s = 0.45 nS, the propagation 
is unstable (upper triangle) except for the initial condition at 
point B2. 



B. Results for a Id Homogeneous Chain of LIF 
Neurons 

Our results for homogeneous chains of LIF neurons 
were qualitatively similar to the above results for chains 
of HH neurons and we do not show the corresponding 
results. Instead of initiating the chain with a current 
pulse, the first LIF neuron of the chain was initiated with 
a synaptic current corresponding to a burst of equally- 
spaced spikes whose times of occurrences were specified 
as initial data (see Eq. (J2J). Similarly, instead of a maxi- 
mum conductance g s being varied, the synaptic coupling 
strength n of Eq. was varied between 1 and 32 which 
we found to span the same kinds of states discussed in 
Figs. Inland [fJl from a 0-spike final state to an initial burst 
that grew without bound in width and in the number of 
spikes. Because the LIF neurons have no intracellular 
calcium dynamics, the number of spikes in the final state 
was primarily determined by the initial spike number. In 
contrast, the intracellular calcium current of a HH neu- 
ron can boost depolarization and maintain propagating 
bursts in cohesion so that chains of HH neurons have 
bigger basins of attraction for stable bursts. 

C. Other Hypotheses and Related Theory 

We discuss briefly here some other hypotheses and re- 
lated theoretical work that might explain the data of 
Hahnloser et al Q and of Mooney 01 but that do not 
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FIG. 6: Plot of the number of spikes observed in a final non- 
transient propagating burst for different initial amplitudes of 
a square current pulse of fixed duration 20 ms and for dif- 
ferent maximum synaptic conductances g s . Different basins 
of attraction are observed, with 0-spike and unstable bursts 
having the largest basins, and roughly equal size basins for 
bursts with 1 to 5 spikes in the final burst. 



involve the propagation of bursts through a feedforward 
network. 

While a feedforward network, especially a synfire 
chain, is one of the simplest concepts that might explain 
sparse precisely-timed bursting, there is much theoreti- 
cal work dating back to the 1980s which shows that a 
recursive network with inhibitory and excitatory connec- 
tions is capable of learning and producing different kinds 
of temporal sequences (32, E3- These networks gener- 
alize the Hopfield attractor model of associative mem- 
ory p3, 43] by allowing non-symmetric couplings between 
pairs of neurons, and by using two kinds of synapses, 
"fast" synapses that stabilize a given network state, and 
"slow" synapses that cause successive transitions between 
the quasistatic network states. Thus there is no difficulty 
in principle for a recursive network to store, generate, 
or learn many temporal sequences of the sort observed 
in HVC. (We note that other classes of recursive mod- 
els are possible, e.g., Huerta and Rabinovich [3i| dis- 
cuss a model with neurons that are strictly inhibitory 
or strictly excitatory, rather than allow any neuron to 
have inhibitory and excitatory connections with other 
neurons.) 

Asymmetric Hopfield-like recursive networks have sev- 
eral attractive features for modeling HVC dynamics. The 
HVC inhibitory neurons [l5j can be incorporated in a 
natural way since temporal sequences are stored via neu- 
ronal connection strengths that typically have pos- 
itive (excitatory) and negative (inhibitory) values. Re- 
cursive networks that generate sequences can evolve via 
a simple Hebbian learning rule 32] from densely intercon- 
nected neurons (see Figs. 2 and 3 of Ref. |15|) and so are 



possibly easier to form during maturation of brain tissue 
than purely feedforward networks. Asymmetric Hopfield 
models are naturally fault-tolerant and error-correcting 
so that details of a stored sequence are weakly affected 
if neurons or synapses are modified, deleted, or added. 
Finally, a single network of this kind is capable of storing 
and generating many different temporal sequences, which 
is consistent with the ability of some songbirds to learn 
and sing many different songs [2j. 

However, further work is needed to determine whether 
existing recursive models (3^ are capable of taking into 
account specific experimental details of HVC such as the 
fact that an HVCra neuron fires briefly just once dur- 
ing a one-second-long motif (a neuron is more likely to 
fire multiple times during a motif if there are recurrent 
pathways), that time interval between successive bursts 
is short (estimated to be about 10 ms [||, this is possibly 
too short for a separation of time scales to exist between 
slow and fast synapses [32^1. that there is a precise align- 
ment of bursts with auditory features during singing 4] 
and during audition [15|. and especially that HVC in- 
hibitory neurons fire tonically while the HVC projection 
neurons fire sparsely during singing 4] and during au- 
dition |15|. (There is no dynamical distinction between 
inhibitory and excitatory synapses in the recursive mod- 
els p.) 

Another alternative explanation for the origin and pre- 
cise alignment of HVCra bursts during singing and dur- 
ing audition is that these bursts are driven by precisely 
timed external inputs. This possibility is suggested by 
the fact that other nuclei such as NIF and Uva are known 
to provide distributed input to HVC, and to all three 
classes of HVC neurons However, some exper- 

iments suggest that this is a less plausible explanation 
than intrinsic generation of bursts within HVC. For ex- 
ample, preventing auditory input to HVC by destroying 
a bird's cochlea (which deafens the bird) or by lesioning 
nucleus NIF does not prevent a bird from singing, nor 
does it cause a bird's song to change substantially over 
time scales of days. Although extracellular recordings 
of an awake singing bird of the kind reported by Halm- 
loser et al Q have not been carried out after deafening or 
lesioning, the fact that the song does not change imme- 
diately suggests that the same pattern of sparse bursts 
should still be observed in HVCra neurons. Since the 
inputs to HVC remain incompletely characterized, ex- 
ternal driving, or some combination of external driving 
and intrinsic HVC circuitry, might be able to explain the 
observed bursts of HVCra neurons during singing. 

A variation of the idea of precise external input driving 
the bursts in HVC is the possibility that part of HVC 
acts as a central pattern generator (CPG) that drives 
the HVCra neurons. For example, the time-ordering of 
bursts in HVCra neurons observed by Hahnloser et al Q 
could arise from a CPG composed of other HVC neu- 
rons such that the CPG connects to the HVCra neurons 
via axons of different lengths and so with different de- 
lay times. (The dynamics of the CPG itself could be 
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explained as a non-symmetric Hopfield model [32j or a 
cyclic synfire chain.) A CPG-based model would differ 
from a synfire chain in that a given burst does not trigger 
causally the next burst like dominoes falling over along 
some path. For example, destroying all HVCra neurons 
that burst at a certain time in a CPG model would not 
stop the firing of HVCra neurons that normally would 
burst a short time later. A recent in vitro slice experi- 
ment 0| suggests that HVC has intrinsic oscillatory dy- 
namics consistent with the existence of a CPG, but it is 
not yet known how these slice data relate to the dynamics 
of HVCra neurons during song. 



IV. CONCLUSIONS 

For an idealized one-dimensional homogeneous feed- 
forward chain of excitatory non-bursting neurons, we 
have shown by numerical calculations that a brief high- 
frequency burst of two to six spikes can propagate in a 
stable way for various choices of parameter values. Previ- 
ous studies of excitatory networks have not addressed the 
propagation of bursts and so have not been directly rel- 
evant for recent experiments concerning the properties 
of nucleus HVC. Previous studies |9( have also mainly 
involved networks of LIF neurons (to reduce the com- 
putational effort) and so have not taken into account 
the more complex dynamics and realistic time scales of 
HH neurons, especially the effects of a calcium current 
which experiments have shown to be present in HVC neu- 
rons [T^ . |37| . Our results show that stable bursts exist 
over a range of parameter values for both LIF and HH 
neuronal models so that the propagation of stable bursts 
in a homogeneous chain is not sensitive to the choice of 
model or of model parameters. 

These results make plausible the hypothesis that the 
sparse precisely-aligned high-frequency brief bursts ob- 
served by Hahnloser et al Q in HVCra neurons during 
singing are intrinsic (do not need external input) and are 
causally connected in that one neuron initiates bursting 
in the next neuron and so on until the end of the chain. 
However, other explanations of the sparse precise bursts 
are possible as discussed in Section IIII CI and further 
experiments, especially using methods that can analyze 
many neurons at once with good time resolution, would 
be valuable in helping to distinguish feedforward chains 
from external driving, recursive networks, or other mech- 
anisms. 

We have not addressed in this paper why HVCra neu- 
rons burst in the first place in, i.e., why information 
is transmitted as bursts rather than as a population of 
spikes. (Sparse firing can aid the learning of sequences as 
pointed out by Fiete et al [43 but this does not require 
bursts specifically.) There is somewhat of a paradox here 
in that the commonly posited purpose of a burst is to 
increase the likelihood of accurate transmission through 
a noisy network [isl |49| , while a similar goal is achieved 
without bursts by using a synfire chain Thus if a 



synfire chain is the correct architecture for HVCra neu- 
rons, it is less clear why bursts are needed, especially 
since HVC has much less noise than cortical neurons 
and HVC synapses are substantially stronger than corti- 
cal synapses [20j . Unpublished calculations by the au- 
thors 0, 0| that extend the present paper to noisy 
heterogeneous synfire chains do suggest that bursts are 
important for achieving robust propagation even for a 
synfire architecture. 

We finish this paper by summarizing some conse- 
quences of the chain hypothesis and of our calculations 
for recent and future songbird experiments. 

1 . A simple consequence of the feedforward nature of 
an excitatory chain is that an experimentalist in 
principle should be able to initiate singing of an 
arbitrary contiguous segment of a motif. For ex- 
ample, if it is possible to stimulate an intermediate 
part of a chain in such a way that a burst begins 
to propagate toward the end of the chain, a song- 
bird might then sing a motif that starts somewhere 
in the middle and that then continues to its end. 
Similarly, physically terminating the chain at some 
intermediate point (say with a local lesion) could 
cause a motif to terminate before its usual end- 
point, or some combination of stimulating in the 
middle of the chain and lesioning later in the chain 
could be implemented in which case the motif could 
start somewhere in the middle and terminate pre- 
maturely. 

Electrical stimulation of HVC via a single extracel- 
lular electrode in an awake behaving bird does not 
initiate singing [5(J, IU H3 but instead resets a mo- 
tif to its beginning (if the stimulus occurs while a 
bird was singing) . The failure to initiate the singing 
of a motif, or part of a motif, does not rule out the 
existence of an excitatory chain but instead could 
imply that a carefully-arranged pattern of input 
spikes might be needed to initiate activity at some 
point in the chain, especially if the network is a syn- 
fire chain for which most neurons in a pool must fire 
in near synchrony for neurons in the next pool to 
fire. Identifying the location of HVCra neurons in 
a chain by, say, optical imaging would be a valu- 
able prior step that could suggest how to work out 
an electrical or photo-uncaging protocol [53 that 
could initiate a chain at an arbitrary point along 
its length, or could suggest how to lesion neurons 
that would terminate propagation of a burst hence 
a motif prematurely. 

2. Our calculations suggest that a homogeneous chain 
based on LIF or HH neuronal models should gen- 
erally have multiple basins of attraction. Different 
kinds of asymptotic stable bursts differing in the 
number of spikes and in their total duration should 
be observed in nucleus HVC for different values of 
synaptic strength and other neuronal parameters. 
Further, hysteresis can be expected such that, for 
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identical experimental conditions, different initial 
states will evolve to different kinds of asymptotic 
bursts. 

Hysteresis under fixed experimental conditions can 
be explored by initiating different kinds of bursts. 
This may be possible if the neurons at the begin- 
ning of the chain can be identified and an appropri- 
ate stimulation protocol worked out (see points 1 
and 3). 

Most nonlinear network models, including asym- 
metric Hopfield models and CPGs, also have mul- 
tiple basins of attraction so this is not a distin- 
guishing property of chain models. However, the 
number and types of basins of attraction, and how 
these basins vary under the influence of pharma- 
ceuticals that alter synaptic strengths, may lead to 
predictions that distinguish feedforward from re- 
cursive models. 

3. Chains of neurons, like the neuronal models that 
make up the chains, are driven dissipative systems 
so that initial states will generally evolve through 
a transient before settling into an asymptotic be- 
havior Q • Properties of a burst such as its num- 
ber of spikes and its set of interspike intervals will 
evolve over time, with the transient time depend- 
ing on the choice of initial condition (how far the 
initial state is from the attractor) and on neuronal 
parameters. Our results such as Fig.0|for an homo- 
geneous HH neuronal model suggest that transient 
times can in fact be rapid, just four or fewer suc- 
cessive neurons. Fig. 2 of Hahnloser et al [j does 
not show evidence of transient behavior since later 
bursts do not seem to be statistically different from 
earlier bursts during the same motif. However, only 
a small number (about 20) of HVCra neurons have 
been sampled to date in singing birds, and neurons 
that burst at the beginning of a motif have not 
been found so further study involving more neu- 
rons would be worthwhile. It may also be the case 
that the songbird brain has evolved in such a way 
as to reduce or eliminate the duration of transients, 
for example by starting a chain with a state close to 
the asymptotic burst, or by using inhibitory neu- 
rons to accelerate convergence to synchronized fir- 
ing as bursts propagate from projection neuron to 
projection neuron |55j. 

4. A one dimensional chain is not robust since propa- 
gation terminates if any neuron in the chain dies. A 
more realistic synfire chain must have at least two 
neurons in a pool to be robust, and the neurons in 
a pool must synchronize their spikes to some ex- 
tent to achieve reliable transmission from pool to 
pool 8]. A synfire chain hypothesis thus predicts 
that there must be more than one HVCra neuron 
bursting at any give time (these are the neurons 
that belong to the same pool) and further that neu- 
rons that burst at about the same time should have 



nearly synchronous bursts. Thus Fig 2 of Ref. Q 
may be incomplete and a high-resolution optical 
study of HVC (or possibly a many-electrode study) 
during the singing of a motif or during the audi- 
tion of a bird's own song [l5l | may reveal multiple 
HVCra neurons firing bursts in synchrony. 

Abeles has pointed out that there is an inverse re- 
lation between the number of neurons needed in 
a synfire pool for reliable propagation and the av- 
erage strength of the synapses between pools ||. 
Using paired intracellular recordings, Mooney and 
Prather [2(j have shown that synapses between 
HVCra neurons are stronger than mammalian cor- 
tical neurons by about an order of magnitude (the 
average EPSP magnitude is about 2 mV compared 
to about 0.2 mV in cortex) so chains in HVC would 
be expected to have about an order of magnitude 
fewer neurons in a pool. Given Abeles's estimate of 
about 50-100 neurons in a cortical pool, perhaps ten 
or fewer neurons might be needed in an HVC syn- 
fire pool. Assuming that the HVC bursts of av- 
erage duration 6 ms are non-overlapping, a zebra 
finch motif of duration 0.6 s would require about 
0.6/0.006 w 100 successive pools to span the mo- 
tif, so a synfire chain in HVC might involve about 
100 X 10 w 1, 000 neurons. This is a much smaller 
number than the 40,000 HVCra neurons estimated 
to exist in HVC. This suggests that a synfire chain 
might be too simplistic an architecture to justify 
the large number of HVCra neurons. 

5. It is worth noting that a synfire chain architec- 
ture can explain how an overall precise timing can 
be maintained in HVC even though there is a 
steady and substantial turnover of neurons in HVC 
throughout the life of the songbird 0, |2^| . Since 
perfect synchronization of a pool is not needed to 
guarantee transmission of information through the 
next pool, a fraction of neurons in a given pool can 
alter their properties or fail without degrading the 
transmission of information or its timing. 

In conclusion, our calculations support the hypoth- 
esis that the sparse bursts observed experimentally in 
HVCra neurons during singing can be understood as the 
propagation of bursts through an excitatory feedforward 
synfire chain. However, further experiments and compu- 
tational studies are needed to confirm this hypothesis and 
to rule out the competing hypotheses discussed in Sec- 
tion [^^21 Especially interesting in the near term would 
be to understand some of the quantitative details of the 
experiments, for example the production of a brief 6 ms 
burst of four spikes corresponding to a frequency of about 
600 Hz that is consistent with the known properties of 
HVC neurons and of the HVC microcircuitry J2jj. Al- 
though our results show that a brief burst can propagate 
in a stable manner, we and others have not been able 
to construct neuronal models, or simple networks of such 
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models, that burst as rapidly and as briefly as actual 
HVCra neurons. 



Note: Toward the completion of this paper, the authors 
received a preprint by Jin et al whose results overlap with 
the present paper (the preprint is related to a recent con- 
ference presentation Q). These authors also investigated 
whether brief high-frequency bursts analogous to those 
observed in HVCra neurons could propagate through a 
purely excitatory network, and they too found that this 
was possible for a range of parameters, thereby support- 
ing the hypothesis that the bursting of HVCra neurons 
could be intrinsic to HVC and causal. 



The main difference of the present paper with the 
preprint is that we explored whether bursts could prop- 
agate under the simplest circumstances of a Id, homo- 
geneous, noiseless, purely-feedforward chain using sim- 
ple neuronal models. Jin et al studied a more com- 
plicated generalized (not purely feedforward), heteroge- 
neous, noisy, purely excitatory synfire chain with many 
neurons per pool. Jin et al also used a more com- 
plicated two-compartment conductance-based neuronal 
model that included a low-threshold potassium channel 
that provided a strong spike-frequency adaptation similar 
to that observed experimentally. (This helped to create 
brief high-frequency bursts.) They also included calcium 
dynamics in the dendrite compartment that could lead 
to the firing of calcium spikes. Jin et al were able to 
show that adding the low-threshold potassium channel 
and calcium dynamics made the propagation of bursts 
through their generalized synfire chain more robust (one 
did not have to tune parameters so carefully to obtain a 
stable propagating burst). They also studied how jitter 
in the burst evolved along the length of the synfire chain, 
which we have also done in a forthcoming paper with a 
somewhat different model [l7j . 



Given that our calculations and those of Jin et al have 
left out the role of inhibitory neurons and were based on 
a limited knowledge of HVC neuronal properties and of 
the HVC microcircuitry (e.g., there is no evidence yet 
for a synfire chain in HVC, nor have the conductances of 
HVC neurons been fully characterized), both calculations 
should be considered as modest but complementary and 
useful steps toward understanding possible mechanisms 
for the dynamics of HVC during singing. 



APPENDIX A: HODGKIN-HUXLEY 
EQUATIONS AND PARAMETERS 

The single-compartment HH model of Eq. had the 
following five representative currents Ii(t): 



-ficak = 9l(v -E l ), 

I-Nn = g-Nam 3 h(v - E Na ), 

Ik = gkn 4 (v - E K ), 

Ic a = gca,rn\hx{v - E Ca ), 

Ikch = 9KC a nl(v - E KCa ,), 



(Al) 



with maximum conductances per unit area y (units of 
mS/mm 2 ) given 

g L = 0.003, .g Na = 1.0, g K = 0.22, g Ca = 0.031, 

(A2) 

and corresponding reversal potentials (units of mV) : 



E L 



-70, £ N a = 50, E K 



-80, 



E Ca 



120, 
(A3) 
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TABLE I: Parameter values for the Hodgkin-Huxley model 
Eq. j2J with currents Eq. lAlfl and gating equations 
Eqs. (EU-lETj- 

The currents Eqs. <|A1|) depend on the membrane volt- 
age v(t) and on gating variables m, h, n, mi, hi, and ri\ 
that determine the probability as a function of time for 
certain channel subunits to be open. All the gating vari- 
ables x(v) except the KCa activation variable ri\ obey 
the evolution equation 



(A4) 



with the asymptotic value Xoo(v) having the sigmoidal 
form 



x 00 (v) = (l + e-(«-«o)M) \ 
The time constant t(v) has the form 

t(v) = t 2 + h (l + e- « v ~ ^/^y 1 . 



(A5) 



(A6) 



for the standard HH channels m, h, and n, and the form 



t(v) = t2 + ti exp 



V - V2 
V3 



(A7) 



EkCs 
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for the Ca channels mi and hi and KCa channel n\. 
Table 1X1 gives the parameter values that we used in these 
various expressions. 

The KCa activation variable n\ obeys the evolution 
equation 



TKCa^- = /([Ca]<) -m, 



(A8) 



where the expression [Ca]i denotes the average intracel- 
lular calcium concentration and the function / is given 
by 



/([Ca],) 



[Ca], 



([Ca], + 0.001)(1 + exp(-(u + 40)/20)) ' 

The intracellular calcium concentration accumulates 
when the low threshold calcium channel is open and cal- 
cium ions flow from extracellular solution into intracellu- 
lar solution. Thephenomenological model for [Ca]^ was 
taken from Rcf. 1501 



[Ca], = ([Ca] i0 - [Ca], - F * (J ca ))/r Co , 



(A10) 



where the equilibrium Calcium concentration [Ca]io had 
the value .00005 mM and the current-density factor F 
had the value 3 mM — cm 2 /mA, and the time con- 
stant Tea had the value 200 ms. The KCa channel was 
chosen to imitate a low-threshold transient calcium cur- 
rent observed in HVC neurons by Kubota [ijj- 



ACKNOWLEDGMENTS 



The authors thank Richard Mooney, Jonathan 
Prather, Melissa Coleman, Michale Fee, Stephen Shea, 
George Augustine, and Paul Tiesinga for helpful discus- 
sions, and Craig Henriquez and Duke University's Center 
for Neural Engineering for providing support for one of 
the authors (M. Li). 



[1] A. J. Doupe and P. K. Kuhl, Annu. Rev. Neurosci. 22, 
567 (1999). 

[2] C. K. Catchpole and P. J. B. Slater, Bird Song: Biological 
themes and variations (Cambridge University Press, New 
York, 1995). 

[3] F. Johnson, K. Soderstrom, and O. Whitney, Behavioural 

Brain Research 131, 57 (2002). 
[4] R. H. R. Hahnloser, A. A. Kozhevnikov, and M. S. Fee, 

Nature 419, 65 (2002). 
[5] A. Reiner, D. J. Perkel, C. V. Mello, and E. D. Jarvis, 

Ann. N. Y. Acad. Sci. 1016, 77 (2004). 
[6] M. S. Fee, A. Kozhevnikov, and R. H. R. Hahnloser, Ann. 

N. Y. Acad. Sci. 1016, 153 (2004). 
[7] D. Jin, F. Ramazanoglu, and H. S. Seung (2005), society 

for Neuroscience Abstracts No. 753.15. 
[8] M. Abeles, Corticonics: Neural Circuits of the Cerebral 

Cortex (Cambridge University Press, Cambridge, Eng- 
land, 1991). 

[9] M. Diesmann, M.-O. Gewaltig, and A. Aertsen, Nature 

402, 529 (1999). 
[10] R. 0§an, R. Curtu, J. Rubin, and B. Ermentrout, J. 

Math. Biol. 48, 243 (2004). 
[11] A. D. Reyes, Nature Neuroscience 6, 593 (2003). 
[12] D. Golomb and Y. Amitai, J. Neurophysiol. 78, 1199 

(1997). 

[13] V. Litvak, H. Sompolinsky, I. Segev, and M. Abeles, Jour- 
nal of Neuroscience 23, 3096 (2003). 

[14] T. P. Vogels and L. E. Abbott, J. of Neuroscience 25, 
10786 (2095). 

[15] R. Mooney, J. Neurosci. 20, 5420 (2000). 

[16] K. Ohki, S. Chung, Y. H. Ch'ng, P. Kara, and R. C. 
Reid, Nature 433, 597 (2005). 

[17] M. Li and H. Greenside (2006), in preparation. 

[18] M. Li, Ph.D. thesis, Duke University (2006). 

[19] M. Kubota and I. Taniguchi, J. Neurophysiol. 80, 914 



(1998). 

[20] R. Mooney and J. F. Prather, Journal of Neuroscience 

25, 1952 (2005). 
[21] P. Dutar, H. M. Vu, and D. J. Perkel, J. Neurophysiol. 

80, 1828 (1998). 
[22] J. M. Wild, M. N. Williams, G. J. Howie, and R. Mooney, 

Journal of Comparative Neurology 483, 76 (2005). 
[23] R. D. Traub, J. G. R. Jefferys, and M. A. Whittington, 

Fast Oscillations in Cortical Circuits (Bradford Books, 

New York, 1999). 
[24] F. Nottebohm, J Neuroscience 22, 624 (2002). 
[25] N. Wang, P. Hurley, T. Pytte, and J. R. Kirn, J. Neurosci. 

22, 10864 (2002). 
[26] Y. Ikegaya, G. Aaron, R. Cossart, D. Aronov, I. Lampl, 

D. Ferster, and R. Yuste, Science 304, 559 (2004). 
[27] A. Leonardo and M. S. Fee, Journal of Neuroscience 25, 

652 (2005). 

[28] K. Doya and T. J. Sejnowski, in Central Auditory Pro- 
cessing and Neural Modeling, edited by Poon and Brugge 
(Plenum Press, New York, 1998). 

[29] H. D. I. Abarbanel, L. Gibb, G. B. Mindlin, and S. Ta- 
lathi, J. Neurophysiol. 92, 96 (2004). 

[30] L. Berger, Master's thesis, Institute of Neuroinformatics, 
University /ETH Zurich (2004). 

[31] P. J. Drew and L. F. Abbott, J. Neurophysiol. 89, 2697 
(2003). 

[32] D. Kleinfeld and H. Sompolinsky, in In Methods in Neu- 
ronal Modeling: From Synapse to Networks, edited by 
C. Koch and I. Segev (MIT Press, Cambridge, MA, 
1989), pp. 195-246. 

[33] R. Huerta and M. Rabinovich, Phys. Rev. Lett. 93, 
238104 (4 pages) (2004). 

[34] P. Dayan and L. F. Abbott, Theoretical Neuroscience: 
Computational and Mathematical Modeling of Neural 
Systems (MIT Press, 2001). 



13 



[35] D. Kincaid and W. Cheney, Numerical Analysis 

(Brooks/Cole, Pacific Grove, California, 1996), 2nd ed. [46 
[36] MathWorks, Matlab, http://www.mathworks.com/products/matl 
[37] M. Kubota and N. Saito, Journal of Physiology 440, 131 [47] 
(1991). 

[38] A. A. Prinz, C. P. Billimoria, and E. Marder, J. Neuro- [48] 
physiol. 90, 3998 (2003). [49] 

[39] J. Huguenard and D. McCormick, J. Neurophysiology 

68, 1373 (1992). [50] 

[40] C. Koch and I. Segev, eds., Multiple Channels and Cal- 
cium Dynamics (MIT, 1998). [51] 

[41] M. Hines, J. W. Moore, and T. Carnevale, Neuron, 

http://neuron.duke.edu/. [52] 

[42] J. Hertz, A. Krogh, and R. G. Palmer, Introduction to 

the Theory of Neural Computation, no. I in Santa Fe In- [53] 
stitute Studies in the Sciences of Complexity (Westview, [54] 
Boulder, CO, 1991). 

[43] J. J. Hopfield, Proc. Nat. Acad. Sci. 79, 2554 (1982). [55] 

[44] M. J. Coleman and R. Mooney, Journal of Neuroscience [56] 
24, 7251 (2004). 

[45] M. J. Rosen and R. Mooney, J. Neurophysiol. 95, 1158 



(2006). 

M. M. Solis and D. J. Perkel, J. Neuroscience 25, 2811 
al2j305). 

I. R. Fiete, R. H. R. Hahnloser, M. S. Fee, and H. S. 
Sueng, J. Neurophysiol. 92, 2274 (2004). 
J. E. Lisman, Trends Neurosci. 20, 38 (1997). 
L. M. de la Prida, N. Stollenwerk, and J. V. Sanchez- 
Andres, Physica D 110, 323 (1997). 

E. T. Vu, M. E. Mazurek, and Y.-C. Kuo, J. Neuroscience 
14, 6924 (1994). 

E. T. Vu, M. F. Schmidt, and M. E. Mazurek, J. Neuro- 
science 18, 9088 (1998). 

D. S. Vicario and H. B. Simpson, J. Neurophys. 73, 2602 
(1995). 

A. Losonczy and J. C. Magee, Neuron 50, 291 (2006). 
P. Berge, Y. Pomeau, and C. Vidal, Order Within Chaos 
(John Wiley & Sons, New York, 1984). 
N. Kopell and B. Ermentrout, PNAS 101, 15482 (2004). 
G. Turrigiano, G. Lemasson, and E. Marder, J Neuro- 
science 15, 3640 (1995). 



